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Abstract 

A recently proposed generating functional allows the construction of the full set 
, of n-point Green functions in QCD at high temperature and at distances larger than 

' 1/gT. One may then learn how the system maintains its thermal equilibrium in the 

quantum field theory approach, i.e. which process compensates for the important dis- 
! sipation due to collisions. This system may be characterized by quantities which have 

I a classical limit. One finds that the fluctuations of the colour field are not gaussian 

CN ■ ones. A comparison is made with Bodeker's classical approach where a random noise 

^ , is the source of fluctuations. 

PACS: ll.lO.Wx, 12.38.Mh ; keywords: quark-gluon plasma, effective theory, dissipa- 
tion 

O ; 

1 Introduction 

Oh! 

^ I QCD at high temperature possesses several energy scales: T, gT, g^T In 1/ (7 ... , where g{T) <^ 
1. Two effective theories have been constructed, each one associated with a specific scale. 
In the HTL effective theory the scale T"^ is integrated out. Its features and conse- 

quences are well established 0, ^. In the ultrasoft effective theory, the scales and {gT)~^ 
are integrated out Those effective theories exhibit semi-classical features and may be 
obtained via different approaches (quantum theory, response function, transport equations, 
distribution functions) [^, |^, ^, At the scale {g'^T)~^, the physics is expected to be 
non-perturbative, as thermal perturbation theory breaks down because of strong nonlinear 
effects. 

As one goes to larger (space) scales, the importance of dissipation at this scale increases, 
as energy is transferred towards smaller scales. In the HTL case, the damping only occurs 
in some kinematical regions (space-like gluons). The quantum theory is mostly constructed 
out of those regions, which may then be reached via analytical continuations in the energy 
variables. In the ultrasoft case, the damping is of the same order as the momentum, and 
it occurs in the whole energy domain. The approaches split into two branches, the classical 
approach ||^, |^, ^, Q and the quantum approach JT], |], which differ in the resulting set of 
classical equations for the effective field. 

The original derivation, by D. Bodeker, is within the framework of the classical HTL 
effective theory. The fields in this theory are split into soft p ~ gT and ultrasoft p ~ g^T, 
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and the soft fields are integrated out to obtain the following set of effective equations for the 
ultrasoft field 

[D,,F^^]{X)=ml [ v''WiX,v) (1) 

{v.D, + C)W{X, v) = v.E(X) + v) (2) 

where v.D^ is a gauge-covariant drift operator with v'^{vo = 1, v), (v^ = !),/„ = / dQ^/An, 
and C is a collision operator local in X, but non-local in v. ^(X, v) is a gaussian white 
noise term, which compensates for the dissipation arising from the collision term, so that 
the system is in thermal equilibrium. The gluons p ~ T take part in the collective motion 
W{X, v), whereas the gluons p ~ gT are exchanged during the collisions that C represents. 

Alternatively, in the quantum approach, one considers the response of the system to 
a weak external perturbation, i.e. to a source j^^'*(X) which induces large-scale inhomo- 
geneities. The response is an induced current and the resulting set of classical 

equations is 

[D„Fn(X)=j7*(X)+j7^(X) (3) 
r;\X)=ml f v^'W'iX,^r) (4) 

{v.D, + C)W'{X,v) = v.E{X) (5) 

where C is the same as in Eq.(H). These equations are obtained from the Schwinger-Dyson 
equation for the two-point function via a gradient expansion. In this approach, the response 
to a short perturbation dies out because of the damping due to C. In order to draw conclu- 
sions about the system in thermal equilibrium, one has to refer to the fluctuation-dissipation 
theorem, which is generally obeyed by the two-point function. 

Morever, this induced current is the generating functional of the n-point response func- 
tions. Recently, a generating functional has been presented |jlO[ that allows one to obtain 



the full set of ra-point amplitudes, i.e. it contains all the information on the (effective) quan- 
tum system at a scale larger than {gT)^^ and smaller than {g'^T)^^ . As the system is in 
thermal equilibrium, there exists a mechanism which compensates for the energy dissipation 
occurring at this scale. In macroscopic physics, the compensation is attributed to thermal 
fluctuations. The purpose of this paper is to characterize, tentatively, the process that main- 
tains the system in equilibrium, and to compare it to the statistical physics approach. 

In quantum field theory at finite temperature T, a doubling of the degrees of freedom 
occurs, compared with field theory at T = 0. There are several formalisms available [Q. 
The Keldysh basis UTTl , a rotated version of the Closed- Time-Path formalism, offers 



two particular advantages for our purpose: (i) the generating functional of ultrasoft QCD 
is simple (ii) there exist n-point Green functions which are strong candidates for a direct 
connection with the correlation functions of statistical physics, at high T. 
Section ^ summarizes several general features of the (unfamiliar) Keldysh basis: (i) the 
relation between one-particle- irreductible (IPI) vertices and connected Green funtions, (ii) 
the relation between Keldysh Green functions and expectation values of ordered product of 
field operators A^{X). A reader familiar with this basis may go directly to Sec. |2.3| , where 
candidates for an identification with correlation functions are selected. Section ^ is devoted 
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to the quantum approach to ultrasoft QCD. At the scale (g'^Tlnl/g)^^ (with \nl/g ^ 1), 
the connected Green functions of interest are deduced from the IPI vertices and these are 
obtained from the generating functional One resuh is that the fluctuations of the field 
Af^{X) at high T are not gaussian. The HTL case is examined as a limiting case, and 
simple explicit expressions are obtained. In Sec.^, those results are compared with those 
from the classical approach, i.e. from Eqs.(|I])(^, where a gaussian noise compensates for 
the dissipation. Section |] contains a summary of the result and discusses an effective action 
for ultrasoft QCD. 



2 The n-point Green functions in Keldysh basis 

In the functional approach to field theory at temperature T, the time contour that enters the 
action has to run forwards and backwards between t = -co and t = +00. When expressed as 
a single time integral, the result is the doubling of the degrees of freedom, as the field ^4+ lives 
on the forward branch and the field A_ on the backward one. In consequence, to an n-point 
amplitude are associated 2" different components. The familiar representation involves 
and A_ and the corresponding Green functions are expectation values of mixed combinations 



of time-ordered and anti-time-ordered products of operator fields. In the Keldysh basis ITT 
the fields are A^ ± y4_, and the Green functions are expectation values of combinations of 
time-ordered commutators and anticommutators of the operator fields. Among them are the 
response functions, i.e. the retarded amplitudes. 

The present state, in the quantum approach to ultrasoft QCD, follows. From the knowledge 
of the retarded n-point functions, consistency condition and gauge symmetry have lead to 
the whole set of IPI vertices in the Keldysh basis. We now need to construct the connected 
Green functions from the IPI vertices, and to identify the relevant ones for our purpose. 



2.1 Relating connected amplitudes to IPI vertices 

It will be useful to obtain the results in two ways: (i) the naive approach is to go from the 
Keldysh basis to the Retarded-Advanced (R/A) basis, where the relation is immediate, (ii) 
the other approach is to extract the prescription from the general Legendre transform in the 
Keldysh basis. 

Our notation is Pi{p^,Pi). In the Keldysh basis, the indices are 1 and 2, the IPI n-point 
vertices are Kij...n{Pi, P2, ■ ■ ■ Pn) and the connected amplitudes are (jJj-...„(Pi, P2, ■ ■ ■ Pn) ■ All 
momenta are incoming X^iLi ~ 0- R/A basis, a leg may be of type R or of type A, 

and the IPI amplitudes are Tra-.-a- In this basis, the relation between amputated and non- 
amputated functions is straightforward. Indeed, to a leg of type R one attaches a retarded 
(dressed) propagator Gr{P), and to a leg of type A an advanced one Ga{P)- 
The connection between IPI n-point vertices in the R/A basis and in the Keldysh basis has 



been studied in Ref. |jT3|, |T0[. Here we just quote the results and make use of them. One 
defines 

iV(P,) = n(pl') + i = icoth^p° (6) 

G,{P) = N{P)[Gr{P) - Ga{P)] (7) 
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where n{p^) is the Bose statistical weight. One relation between the bases is 

i^21...l(Pl, P2---Pn)= TRA-AiPl, P2 • ■ • Pn) (8) 

SO that 

G 21- ■ -liPl ^ P2 ■■■ Pn) = Gr{Pl)Ga{P2) ■ ■ ■ Ga{Pn) ^RA-a{Pi, P2- ■ ■ Pn) 
G^21...l(Pl,P2---Pn) = GriPl)GaiP2)---GaiPn) K21...liPl,P2---Pn) (9) 

i.e. the relation is now written in the Keldysh basis. Similarly 

i^221...l(Pl, P2 ■ ■ ■ Pn) = TRRA...AiPl, P2 ■ ■ ■ Pn) 

-L ARA---A J- RAA - A (10) 

G^221-- l = GaiPs) ■ ■ ■Ga{Pn)[ Gr{Pl)Gr{P2)^ RRA --A 
-N{P{)Ga{Pl)Gr{P2) -L ARA --A - N{P2)Gr{Pl)Ga{P2) -L RAA---A ] (11) 

^221--1 ~ GaiP-i) ■ ■ ■Ga{Pn)[ C^. (Pi ) Gr (-P2)-^221--l 

+Ge(Pl)a(P2) Kui-1 + Gr{Pl)G,{P2) K2U...I ] (12) 

One notices the occurrence of Gc{Pi) , which will play an important role later on. Similarly, 
the next relation is 

^2221...! = Ga{P4) " ' ■G„(P„) { (A )a (^2)^, (^3)^^2221...! 
+ [ G,iP,)GriP2)GriP3) Ku21-1 + (1 ^ 2) + (1 ^ 3) ] 

+ [ a(Pi)G,(P2)a(P3) 7^211...! + (1 ^ 2) + (1 ^ 3) ] } (13) 

where (1 ^ 3) means the term obtained in the exchange of all variables of legs 1 and 3. The 
denominators of the propagators Gr{P) and Ga{P) are respectively (po + ie)^ — — Ilr{P) 

and (po - le)^ - - n„(p), with n,(p) = n;(p). 

Alternatively, one may obtain those relations via the Legendre transform in the Keldysh 
basis (see Sec .12. 21). The result is, for n > 3 



G'r^k-n = Mu'{Pl)Mjy{P2) ■ ■ ■ M^r^'iPn) K.^.^'-n' (14) 

where the 2x2 matrix M{P) is 



M(P) 



GaiP) 
a(P) Gr{P) 



(15) 



The zero of the matrix M ensures that -ft'ii.-.i = leads to G11...1 = 0. Both are important 
features of the Keldysh basis. Also, the retarded amplitudes G211...1 are simply related to 
K211...1. The 2-point function is written as a 2x2 matrix 

aiG{P)aiH{P) = 1 , M(P) = G(P)ai (16) 



4 



H21 = (po + -p^~ n,(P) , Hu = {po - ^e)^ - - n,(P) 

Hn = , H22 = NiP)[H2i-Hu] (17) 



(Ti is the Pauli matrix, and it arises from the fact that, in the Keldysh basis, the source term 
in the Legendre transform reads (see Eq.(^)) 

J J,A2 + J2A, = J {A, A2) ^1 ( ) (18) 
i.e. the generating functional Ty is a function of aiJ. From Eq.(|16]) one deduces 

G2I = Gf = , G12 = Ga = 1 G22 = Gc = —GrH22Ga (19) 

so that Eq.(|3) translates into Gc{P) = N{P){Gr - Ga). 

For n > 4, it is well known that the connected Green functions contains other terms in 
addition to those written in Eqs.(p!^) ([T3|). They arise from one-particle-reducible pieces, i.e. 
from IPI vertices linked with dressed propagators. For example, for the 4-point case, two 
3-point vertices are linked with one propagator, and the prescription is 

Gf^i = Mu>{Pi)Mjj.{P2)K,j.„,{H~X,„,{Pi + P2)Km'kn'Mk>k{P^)Mvi{Pi) 

+ (2 ^ 3) + (2 ^ 4) (20) 



2.2 The Green functions in position space 

As the Keldysh basis is unfamiliar, since rarely used, we summarize the properties that are 
needed to fully understand the physical features of our results on ultrasoft QCD. Those 
properties are taken from Ref.ll^] (our definition of Gc{P) differs by a factor of 2, see ||13|| ). 



The same generating functional can generate n-point functions in different representations 
of the Closed-Time-Path formalism, provided that the external source term is expressed in 
the corresponding functional arguments. This source term writes 

/[j+(X)A+(X)-j_(X)A_(X)]= / + (21) 

with Jx = dt J d^x, where A+ is the component of the A field along the forward branch 
of the time contour. A- the component along the backward branch, and 

Ai = (A+ + A_)/2 , A2 = A+-A^ ,Ji = (j++j_)/2 , J2=J+-J- (22) 

In the Keldysh basis, the Legendre transform reads 

W{JuJ2)= [ [J1A2 + J2A,] + T{Ai,A2) (23) 

8r sr ^ dw , dw , 

smx) = ■ arr wr ' wr^' 

00 1 „ 

r(A,A2) = E - / i^^,7...n(Xi,---X„) A(Xi)A,(X2)---A„(X„) (25) 



5 



dW 

Since A'li...! = 0, T{Ai,A2) (and all its derivative with respect to Ai only) vanishes for 
A2 = 0. Setting J2 = amounts to A2 = 0. Two consequences are important features of the 
Keldysh basis: (i) W{Ji, J2 = 0) = , (ii) = is A+{X) = A_{X), i.e. after differentia- 
tion of W, setting J2 = says that the field A takes the same value on both branches of the 
complex time path. 

The comparison of the Green functions in the bases {A^,A_) and (^1,^2), lead to the 



following results The amplitudes G211 - i are the usual retarded Green functions, which 
involve commutators of the field Af^[X). A functional derivative with respect to Ji generates 
a commutator, a derivative with respect to J2 generates an anticommutator. For example 

G2i = -2^(1,2)([A(1),A(2)]) , G22 = -^({A(1),A(2)}) (27) 
G211 = Hmi,2,3){ [ [A{1),A{2)],A{3)] ) + (2 ^ 3) 

G221 = {e{W){ [{A{1),A{2)},A{3)] ) + 0(1,3, 2)( [{A(l), A(3)}, A(2)] )) + (!- 2) 

G222 = H)'\eil,2,3){ {{A(l), A(2)}, A(3)} ) + 5 perm (28) 

where 0{i,j, k) = 6{i,j)6{j, k), and one sums over all permutations of 1, 2, 3 in G222 • More 
generally 

G22-2= E ^-^^^(l,2---n))({{A(l),A(2)},---A(n)}) (29) 

perm 

2.3 The high temperature Umit 

One considers our case of interest, i.e. an effective quantum field theory for ultrasoft (or 
soft) gluons where po ^ T, then the Bose factor 

N{P) = n{po) + 1 = 1 coth ~ - » 1 (30) 
2 2 2 Pq 

The occupation number J\f is so large than A/" + 1 ~ A/", and one expects the commutators 



of the fields to be negligible in front of anticommutators, in the classical limit [TT[]. One 
concludes that the amplitudes G22-2 of Eq. (PP]) should correspond to the correlation functions 
of the statistical classical theory, i.e. the connected amplitudes G22...2 should correspond to 
the cumulants of the correlation functions and they should characterize the fluctuations. 
The overwhelming role of the amplitudes with only indices 2 is explicit in the comparison 
of the factors that multiply K21...1 in Eqs.(|9D (p!2D (|T3|). Indeed, each extra index 2 brings an 
extrafactor T/po, since this index appears with a factor Gc{P) ~ {Gr{P) — Ga{P))T/po. 

In Appendix 0, we discuss the generalization of the fluctuation-dissipation theorem to 
the n-point functions' case. 
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3 Ultrasoft QCD 



As summarized in the introduction, ultrasoft QCD is an effective theory for gluon momenta 
p -C gT <^ T. The gluons p ^ T take part in a collective motion, whereas the gluons p ~ gT 
are exchanged during collisions which change the direction v of the collective motion. This 
section is devoted to the quantum approach which uses the response of the system to an 
external perturbation (see Eqs.(D(|) and comments that follow). At a space scale larger 
than {gT)~^, the dissipation is important, and our aim is to characterize the process that 
maintains the system in thermal equilibrium in this approach. We will examine the IPI ver- 
tices in the Keldysh basis and the connected Green functions G22 2 of ^^is quantum theory, 
either at the scale p ~ g^T In 1/g (i.e. an intermediate scale between gT and g'^T), or at the 
scale p ~ gT. 

The equilibrium condition sets constraints on the n-point gluon amplitudes. The ultrasoft 
effective theory satisfies them in an economical way in the Keldysh basis: (i) all IPI ampli- 
tudes with more than two indices 2 vanish, (ii) there is a generating functional jfj,{X, A) for 
the amplitudes with one index 2 [0], and another functional JC^^{X, Y, A) for the amplitudes 



with two indices 2 ||l0[. = is the response of the system to an external perturbation, 
and JC^i, would vanish if there were no dissipation. They are both expressed in terms of one 
Green function 

t {v.D, + C) GretiX, Y- A- V, v') = 5(4) (X - Y) 5s, (v - v') (31) 

C Gret{X,Y■,A■,^r,V)= f C(v, v")aet(X, F; A; v", v') (32) 

with GretiX,Y) = for Xo < Fo , = JdnjAir , v^'ivo = l,v),(v2 = 1) , and 
is the covariant derivative Z)^ = + i(yfT'^A^j(X) with T"" in the adjoint representation. 
The collision operator C is symmetric in v and v' and has positive eigenvalues of order 
^f^Tln 1/g, except for a zero mode. Its precise form is irrelevant for our purposes. To leading 
log accuracy, the form of C has been derived in the different approaches |, H, |^, H, H. 



Attempts to go beyond this may be found in Ref. |T5|, |T6|. 

Eg . (|3l|) for Gret may be solved by recurrence as an expansion in powers of C(v, v'). It 
leads to the following interpretation. A collision changes the direction of v. Between two 
collisions, Gret{X,Y; A;\,v') describes a straightline propagation (in direction v) in the 
background non-abelian field [|T^. Alternatively, the solution to Eq.(^) may be written 



as an expansion in powers of A^, i.e. between two interactions with the A field there may 
occur many collisions which change the direction of v. 

r^iX;A) = ml[-g,^A^,iX) 

+1 ! j d'Y V, Gtt{X, Y- A- V, v') dy,Xv' .A\Y)) ] (33) 



JCZ{X,Y;A) = miT 



f [ v,Gtt{X, Y- A- V, v'X + vM'rltiX, X; A; v, v')< ] (34) 
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where m|, = Functional derivation gives the n-point IPI vertices 

Qn-l 

,»-A-.,.,(.Y; .Y. . . . AVO = UX: A) U., (35) 

K22{X, Y) = /C^j,(X, Y] A = 0), K221 is obtained from one functional derivative of /C^^jy and 



so forth. An important identity is [|1J, [16 

dA'dXi) 

= I aei(X,Xi;A;v,v")^7rV',aet(Xi,F;A;v",v') (36) 

The bare vertices should be added to those effective vertices. However, they are subleading 
for momenta p <C gT, as it will be shown explicitely in Sec.^. In the following, we concentrate 
on these IPI vertices. Such an expansion of Gret{X,Y\A\\,w') in powers of A is only 
meaningful for distances such that {v .D C)^^ can be expanded in powers of A. We need to 
assume In l/^f 3> 1 so that the scales g'^T In l/(yf and g^T are well separated. For an ultrasoft 
field gA ~ g^T, and we shall sit at the scale dx ~ g^Thi 1/ g ^ gA. As an aside comment, 
these IPI vertices may be used to take into account the contribution of loops with momenta 
g^Tln 1/(7 to the amplitudes at the scale g^T. The resulting term in the collision operator 
C is local in X only for In l/gf ^ 1. 

3.1 The ultrasoft case 

In this section, expressions for the connected Green functions are derived for momenta p ~ 
(yf^Tln 1/(7, mostly concentrating on those with only indices 2. 

For the 2-point functions, the functional derivatives lead to the following results 

7^21 (P, -P) = K2i{P) = ml6-' [ f v,{v.P + iCyX (po) - g^oQuo ] (37) 

Ki2 is obtained from K21 with the following substitution fi P —P. 

K22{P) =mlT S'^' f v^iv.P + iCyhl + v,{-v.P + iC)-\ (38) 

J v.v' 



so that K22 = {K21 — Ki2)T/pq, which is relation ([T7|) in the approximation N{pq) ^ T/p^. 
K22 exists in the whole kinematical domain (po;P)) with the constraint PqiP <^ gT. Ward 



identities relate the components |T0[ P^Kl^i = = Pf^K 22 ■ explicit expression for K21 



exists in terms of continued fractions |T5| , |T6|. 



The definition of the retarded and advanced propagators Gr{P) and Ga{P) requires the 
choice of a gauge and the projectors V^^{P) on transverse and longitudinal components. 
They are the standard ones and they will not be needed. The denominators are respectively 
P2 - n,(P) and P2 _ n„(P) with K2i{P) = n^(P), Ki2{P) = n^(P). The following useful 
relation is gauge-independent 

6-22 = = - [ G'r(p) - Gr(p) ] 

Po 

= ^ Gr(P)[n.p.(P) - n,,.(P) ]Gr(P) (39) 

Po 
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n,(P) = n;(P) = n,(-p) 



(40) 



Turning to the 3-point amplitudes, all of them may be expressed in terms of the same 
quantities T{i,j,k). If (Pi,/i, a) refers to the momentum, polarisation, colour of gluon 1, 
(P2, u, h) refers to gluon 2, (P3, p, c) to gluon 3, one defines the "time-ordered" quantity 

r(l, 3, 2) = ml [ v^iv.Pi + tC)-\'^ {-V.P2 + 2(7)- V (41) 

T(1,3,2)* = T(2,3,1) (42) 

Pi + P2 + P3 = 0. The order in v space and in colour space {if"''^^) follow the "time-order". 
The IPI amplitudes have either one index 2 (they are the retarded amplitudes), or two 
indices 2, as ^^^222 = 0. For example, 

i^2ii(Pi, P2, P3) = tr'^ni, 2, 3)(-p°) + tr'm, 3, 2)(-p°) (43) 

K22i{Pi, P2, P3) = ^^'T [ r(l, 3, 2) - T(2, 3, 1) ] (44) 

The minus sign in K221 comes from the colour factor — — f"''^^, and there is a cancellation 
between the real parts of the T integrands, i.e. -R'221 is a real quantity, symmetric in the 
exchange of all variables of gluons 1 and 2 (momentum, Lorentz, colour). One sees explicitely 
that the existence of the amplitudes with two indices 2 is related to the presence of the 
damping caused by C. For C -C Pi, K221 is proportional to if°''^'^{iC). Ward identities 
connect these functions to the 2-point ones, in particular (the colour is factored out) 



zP,,K^^r{P,,P2,Ps) = K'2^{P2) 

-tPs,K^27iPi,P2,P3) = K^2{Pl)-K^2iP2) (45) 



An example in position space is the case fi = u = p = 

K^l{X,Y,Z) = -ir'T I f^f^ r°o°(Pi,P3,P2) 
221V ; J J (27r)4(27r)4 ^ ' ^' ' 

sin(Pi + P2).{-Z + (X + Y)/2) sin(Pi - P2).(X - Y)/2 (46) 

antisymmetric in the exchange X ^ Y 

One now deduces from Eqs.(^3D (^4]) the expressions for the connected Green functions. 
The needed relations were written in Sec.|2.1|, i.e. Eqs.(|^)(|T2|)(|TBp, together with the approx- 
imation N{Pi) ^ T/p^. One obtains 

G^211 = GriPl)GaiP2)GaiP3) 1^211 (47) 



G^221 = GriPl)Ga{P2)Ga{P3) T(l, 3, 2) T ^r^" 

+G',(Pi)a(P2)G',(P3) r(2,3, 1) T tf''^ 
+G,(Pi)G,(P2)G,(P3) T(l,2,3) i-pl) ir'^ 

+Ge(Pi)a(P2)G,(P3) T(2, 1, 3) {-pD if'^' (48) 
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One notices that the following product always appear 



Compared to K221, i.e. Eq.(^3|), G221 has two extra T terms, each one multiplied by a factor 
Gc{P). All terms are of the same order as Gc{P) carries a factor T/pq. Endly, from Eq.(|l^) 

^222 = [ Gr{Pl)G,{P2)G,{P3)K2n + (1 ^ 2) + (1 ^ 3) ] 

+ [ ^(^1)^,(^2)^(^3)^^221 + (3 ^ 2) + (3 ^ 1) ] (49) 

G'222 = T G^iPs) [ G,(Pi)a,(P2)r(l,3,2) 

-G„(Pi)a(P2)T(2, 3, 1) ] + (3 ^ 1) + (3 ^ 2) (50) 

where (3 1) means a term obtained by the exchange of all variables of gluons 3 and 1 (the 
colour order follows the time-order). One goes from Eq.(^) to Eq. (pOD by noting that the 
quantity T(l, 3, 2) occurs in 7^211 and in i^22i so that 

- pI Gc{P2) + T Gr{P2) = T Ga{P2) (51) 

Alternatively, 6*222 iiiay be written in terms of the retarded functions (see Appendix |^). 
One sees that while the IPI amplitude -ft'222 = 0, the amplitude G222 differs from zero. It 
is completely symmetric in the exchange of all variables of any pair of gluons. For C <^ Pi, 
G222 is proportional to if°'^'^{iC)^. 

One may convince oneself that the antisymmetry of the colour coefficient f"''"^ does not 
force the symmetric G222("^'^'^) ^'^ vanish. A simpler symmetric example is the sum 
(7^221 + ^212 + ^122) (AT, Y, Z) for the case ^ = u = p = 0, with K^^^ given by Eq.(||). This 
combination is antisymmetric in the exchange of any two space-time variables, it vanishes 
when two point coincide, or at mid-points such as Z = {X + Y)/2, however it does not vanish 
in general. As it has been argued in Sec.p73|, 6*222 (AT, F, 2') only involves anticommutators 
of the field A^, and it should correspond to the 3-point correlation function of statistical 
physics. Then, G222 ("^5 ^5 -^) is the third cumulant of these correlation functions. As it 
does not vanish, one concludes that the fiuctuations of the Afj_ field do not have a gaussian 
character. 

One now considers the 4-point function 6*2222- The terms that come from the 4-point IPI 
vertex have a structure similar to the 3-point one, i.e. Eq . (pOD , through the same process 

m 

G\ill = T G,{P2)G,{P^) [ a(Pi)G',(P4)r(l,2,3,4) 

+Ga{Pi)Gr{PA)T{A, 3, 2, 1) +(2^3)] + 5 terms (52) 

where -|- 5 terms means the five other ways of choosing a pair among (1, 2, 3, 4), and 

r(l,2,3,4) = f^bmjmcd I 

J all V 

v^iv.Pi + tC)-%iv.{P^ + P2) + tCy'v;i-v.P, + tC)-'v': (53) 
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Ward identities relate T(l,2,3,4) to T{i,j,k). For C ^ Pi, 6*2222 is proportional to 
(iC)^- The contribution G'2222' from a pair of 3-point vertices joined by one dressed prop- 
agator is obtained from Eg . (^Of) . The resulting expression is lengthy, there are terms with 
GriP^ + Pj),GaiPi + Pj) and GciPi + Pj). For C < Pi, it is proportional to (iCf. G'2222 
's structure is different from G'2222 's one and the sum does not vanish (they cooperate in 
current conservation in some circumstances). One concludes that the cumulant of order 4 of 
the correlation functions does not vanish. 

For the n-point case, the term that comes from the n-point IPI vertex is an immediate gen- 
eralization of the 3- and 4-point cases, i.e. Eqs.(pOD (^). To conclude, the ultrasoft gluons 
exhibit an infinite series of correlation functions' cumulants G22 2- These cumulants char- 
acterize the fluctuations that maintain the effective quantum system in equilibrium, at an 
intermediate scale between {gT)~^ and (g'^T)^^. 



3.2 The HTL case 

The HTL theory is the effective theory for the soft modes p <^ T , i.e. when the scale T^^ 
is integrated out. Here, the damping only occurs in the kinematical domain of space-like 
gluons. This case is a limiting case of the above formalism, as it corresponds to p ~ gT ^ 
C ~ g^T In 1/(7 , so that C may be approximated as 

C(v,v0^e55.(v-v') , e>0 (54) 

(v.P + iCy^ ^ P.P.— -iTT 6(v.P) (55) 
v.P 

This case allows more explicit expressions for the amplitudes related to dissipation. For the 
2-point function, the result is well known as Landau damping 0]. From Eq. (p7|) with (|55|), 
one obtains for P'^ = — > 0, for example 

K^iiP) = mU^fp, Lipo,p)=ml{^rp, ^\n^ (56) 
p p 2p po + p 



and from Eq. (|38|) with 

T 
Po 

For example 



^22' = -(Hr - n,) = mlT f (-2m) 6{v.P) v^v, (57) 
Vn Jv 



= -m^zTT^ e(-n , K^^ = ^-^Kll (58) 
p p"^ 

consistent with P^K22 = 0- Moreover, Gc = {Gr — Ga)T jp^ is, up to a factor, the transverse 
or longitudinal spectral function of the soft gluon H]. 



For the 3-point function, T^^p{1, 3, 2) is given by Eq.(^) with the approximations (p^) (p5|) 
It has been expli cited in Ref.[0. For the case n = u = p = and for > 0, P2 > 

-Wr»"(1.3.2)^M(P„P.)^^,„^^ (59) 
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A = (P1.P2)' - PfPl = A(Pi, P2) (60) 

M(Pi, P2) has branch points at = 0, Pa^ = and A = 0. At the branch point Pf = 0, M 
gets an imaginary part ^in 

For space-hke indices /i, i/, p, the tensor structure is decomposed on the vectors n = pi x 
P2 , n X pi , n X P2 and one obtains the result that for Vi, or ViVj, or ViVjVk, the integrals 
are linear functions of M(Pi,P2) and of L(Pi),L(P2) (defined in Eg. (|56|) ) with rational 
coefficients [|]. 

7^221 carries a definite prescription for the imaginary parts. From Eq.(^lD with (^) 

ir2°2°? = mlnr'^J^=[ e{-Pl) - e{-Pi) ] (61) 



A(Pl,P2) 

with A defined in Eg . (|60D . For K22i{Pi, P2, P3) to differ from zero, one needs either > 
0, P| < or P2 < 0, P| > 0. 



^221 = A[ P2(n X pi). - P?(n X P2). ]K°2°? 

+m^T^[ (n X p,),-e{-P^) + (n X p2),-^(-P|) ] (62) 



consistent with 



-^[p',K^T,-p\K^^i] = K^^iP.). 

Turning to 6*222 5 the term explicitely written in Eq.(|50D differs from zero if P| < and 
Pf < (or P2 < 0). For example, for the case /i = z/ = p = 0, 6*222 differs from zero in 
narrow domains of the plane Pi + P2 + P3 = 0, such as the domain < Pi, IP2I > P2 with 
bsl < (p?+pi + 2pi.p2)i/2. 

A similar treatment may be made for the 4-point function. T(l, 2, 3, 4) of Eg . (^31) may be 
expressed in terms of M(Pj, Pj) and L{Pk) functions [Q. Its imaginary part may be obtained 
from Eg.(|53D with Eg. (^5]) , and one sees that it is nonvanishing for Pf < 0, or P| < 0, or 
(Pi + P2)^ < 0. One concludes that the soft gluons exhibit the same type of correlations as 
the ultrasoft gluons. 



4 Classical correlation functions in ultrasoft QCD 

In this section, a comparison is made, at the scale p ~ g^Tlnl/g, between the Green 
functions G22 2 '^^ the guantum approach and the corresponding guantities of the classical 
approach to ultrasoft QCD. D.Bodeker's method was to integrate out the scale {gT)~^ in 
the kinetic eguations associated with the field modes p ^T. The field modes gT < p < T 
were split from the softer modes, their eguation of motion was solved and plugged into the 
eguation for the softer modes. A thermal average was performed over the initial conditions 
which enter that solution. The resulting classical eguations for the ultrasoft field modes 
p <^ gT are expressed in terms of an auxiliary colour field W{X, v) 

(v.D., + C)WiX, v) = v.E(X) + e(X, v) (63) 
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[D^,FniX)=ml [ t;^iy(X,v) (64) 

where ^(X, v) is a gaussian white noise with vanishing average, and with correlation function 

2T 

(( UX, v) UX', v') )) = —Civ, v') 5,, 5\X - X') (65) 

where (( )) means an average over the random noise. The noise arises from the average over 
initial conditions for the modes gT < p < T that are integrated out. This noise maintains 
the ultrasoft modes in thermal equilibrium, as it compensates for the dissipation due to C. 

A related idea has been presented for the HTL theory, i.e. for the field modes p <^ T 
r?! , |18| . The classical equations for this field are Eqs. (^) (|6^ ) with (7 = and ^ (X, v) = 0. 



A summary of the results are: (i) in the abelian case (or the abelian approximation to QCD) 
the average over the initial conditions of the W field produces a fluctuating piece C(-^) iii 
the equation for the Ai field (in the gauge = 0), (ii) the average over the W field's initial 
conditions amounts to a gaussian functional integral, (iii) the fluctuating C{X) generates the 
Landau damping part of the Ai field's 2-point function. 



In this section, we solve the set of equations (^)(0) as an expansion in S,-, and we compare 
the noise average of the 2-point and 3-point functions of the A^ field with the correlation 
functions G22 and G222 of the quantum approach (see Sec.p.3|). 

Eq. (|63|) is solved by means of the Green function defined in Eq.(^). With initial condition 

W = Q ait = -00, 

W (X, v)=i( I GretiX, Y- A- V, v') [ v'.E(r) + ^(r, v') ] (66) 

JY Jv' 

with Jy = J d^Y and J^, = J dQ^/fAii. Then, Eq.(^) for the ultrasoft field turns into 

[D„ F^-^UX) = j:(X, A) + j-,(X, A) (67) 
where ja{X, A) may be written as in Eq.(^) and 

jl^{X,A)=mlt f f Gtt{X,Y-A-v,V)i{Yy) (68) 

JY Jv,v' 

Eq.(|67|) may be solved by iteration as an expansion in ^, at the scale X ~ [g'^T In 1/g)^^. 
First, one expands the currrents j'^{X,A) and j^^^{X,A) in powers of the gauge field. This 
expansion is valid, i.e. dx ~ g^Tlnl/g ^ gA ~ g'^T provided Inl/g ^ 1. From Eg . (|35|) , 
formally 

j^(X, A) = W A, + I T'^-'P A,Ap + --- (69) 
and a similar one for j|^(X, A). One, then, writes Eg . (|67|) as 

{d^g^"" - df^d" - W") A, = j[(X, A) - gF^'A, + ^T^^'PA^Ap + ■■■ (70) 

To lowest order, one sets A = in the right-hand-side of Eg.(^). The retarded dressed 
propagator G^^^X, Y) is an inverse of {d'^g^.u — d^'d'^ + Up-'-') , and the lowest order solution 
to Eg.(|7Dp is (with initial condition A = at t = —00) 

AW^(X) = G'r'(X,Xi) mlt vi^, a.(Xi,X2; vi,V2) ea(X2, V2) (71) 
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with the convention that a repeated space-time, or v, variable means a summation over this 
variable, and with the definition 

GriX,, X2; vi, V2) = GretiXi, X^; A = 0; vi, V2) (72) 

One now considers the average over the random noise. (( A^^'^{X) )) = from {{ ^ )) = 0. 
The computation of (( A^^'>'^{X)A^^^'^{Y) )) is instructive, and the steps are detailed in 
Appendix 0. The result is 

(( A^^Hx)AP''iY) )) = / ^ e-^-(^-^)5., .Gr(PO (73) 

with Gc{P) defined in Eg. (|39|) . Since Gc{P) = Gc{—P), this 2-point function is symmetric 
in the exchange X ^ Y. 

A result equivalent to (^) was obtained in the original paper on ultrasoft in the static 
case, and also for the HTL case |jl8|, 0. One may compare the expression ([73|) with the 
quantum theory result, Eqs.(P7|)(P^), and one sees that 



A(^^^{X)A'^>{Y) )) = G^,,{X,Y) = -'-{ {4°)^(X) , Ar^(F)} ) (74) 

One finds an exact correspondence between the (noise-averaged) classical correlation func- 
tion and the matrix element of the field operators' anticommutator. All the higher order 
cumulants of A^^'^{X) vanish, since, from Eq. (|7lD , A^^'^{X) is linear in the gaussian noise 

Turning to the next iterated solution to Eq.(|7D, A^°)^(X) + A^^)''(X), it is obtained with 
the substitution of A^^^'^{X) to A^{X) in the following three terms of the r.h.s. of Eg . ([70|) : 
the terms involving the bare 3-point vertex, the effective vertex T'^"'^, and the second term 
in the expansion oi ^{X, A). These terms are, then, bilinear in ^ and such that 

(( AW'^I^) )) = since UJbc = (75) 
{{A'^:HX)Af'>^{Y))) = since (( ^O) = (76) 

The next step is to compute (( A^^''^ (X) A^^^'^ (Y) A^^^f (Z) )) and to compare the symmetric 
quantity 

Gs = {{ A«^(X)4°)''(r)4°)''(Z) )) + (X, fi,a^Y, u, h) + (X, /x, a - Z, p, c) (77) 
with the quantum result 6*222- O'^^ defines the Fourier transform 

{2nY6\P, + P2 + P3) AZ'iPu P2, P3) (7J 



(( A«-(X)Ar(F)Ar-(^) )) = / n ^ e-(^^-^+^^-^+^3.z) 



One considers successively the contribution of the three terms. With the bare vertex contri- 
bution 

^(i,«)M(X) = -Gr'(X,XO gU, Af^'^iX,) [ 9,, 4°)^ - S.Ajf (X) ](Xi) (79) 
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one obtains (see Eq. ( |105| ) in Appendix ^ 



^r^"^(^i,^2,P3)=^7Gr'(Pi){ 

^fabc Gr{P2) [ Ps.G^ ,.'(^3) - Pa^'G^ .(A) ] + (^2, h ^ P,, p, c) } (80) 
With the effective vertex contribution 

= Gr'iX,X,) |rS,^.(Xi,X2,X3) Ar''(^2)Af''(X3) (81) 

the result is 

■4r^^^ = -G'r'(Pl){Gr'(P2)Gf'(P3)frS,,,(Pl,P2,P3) + (P2,i^,&-^P3,P,C)} (82) 

with the correspondence between the notation in the R/A basis and in the Keldysh bases 
ri?,AA(Pi,P2,P3) = K2ii{PuP2,Pz), and with the exphcit expression for i^aii- The 

symmetry of -R'211 in the exchange of the variables 2 and 3 leads to an identical term in 
Eq.dD so that 2x g/2 = g. 



Comparing the contributions (a) and (/5), i.e. Eqs. (pO|) and (|82|), one sees that the bare ver- 
tex' one is subleading for momenta p ~ g^T In 1/(7, as the effective vertex carries an extra fac- 
tor of order m'jj/{v.P)'^. One now compares Eq.(P^ with the form ( ^5D for the quantum corre- 
lation function 6*222 ' sees that A^^\Pi, P2, P3) is the term G r (Pi) Gc{P2)Gc{P3)K 211, 
so that the symmetric combination defined in Eq. ( [77| ) corresponds to the three terms, in the 
form (|49|), that have one Keldysh index 2. (In Sec.^, g has been factored out of the vertex, 
see Eq.(|5D) 

Endly, the term arising from the second term in the expansion of j^{X, A) is (see Eq. (pB]) ) 
A^^)'^(X) = Gr'(X,XO ml ig i^, ^^^'^'(X^) 

X Wi^/^(Xi,X2, Vi, V2)t;2;.'^(^2,^3, V2, V3) ^g(X3, V3) (83) 

After some computation, one obtains (see Appendix ^ and Eq. ( |10^ )) 

•AZc^'^ = -Gf{Pi){ Gf{P2)Gf{P,) g tUc [ r(l, 2, 3)+S]+ (P2, u, b ^ P3, p, c) } (84) 

with T{i,j, k) defined in Eq.(^), and where S" is a term that will disappear in the symmet- 
ric combination ([77|) (see Appendix 0). One, then, sees that the symmetric combination of 
A^"'^ terms builds up the six T{i,j,k) terms, as they are present in the sum of the terms 



with two Keldysh indices 2 of the form (|49|). This comparison exhibits the correspondence 
between the terms arising from the noise in the classical equation and the terms related to 
the dissipation in the quantum approach. 



To conclude, the symmetric combination ([77|) of noise-averaged approximate classical 



solutions lead to the same value as the matrix element of the combination of anticommutators 
defined in Eq. (P^) , as suggested in Sec. |2.3| . This result brings comfort as, somehow, it was 
expected. Indeed, from the quantum theory point of view, the lack of freedom at the level 
of the 3-point functions is the same as for the 2-point ones. This is best seen in the R/A 
basis where the set of 3-point functions is made of the three retarded ones, such as Traa, 



15 



and their complex conjugate, the advanced ones. As a result, G222 is, in fact, determined by 
the response functions, just as 6*22 is (for the general relation, see Appendix 0). 

We have not pushed any further the comparison between the classical and quantum 
results, due to the proliferation of terms in the classical approach. One may trace the terms 
that are present in G2222'- The third term in the expansion of j^{X,A) will lead, in the 
4-point correlation function, to the terms of the type -R'21115 and the terms of the type -ft'2211 
appear in the same expansion of j^{X, A). 

5 Conclusion 

We have followed a quantum field theory approach to ultrasoft QCD. From the knowledge 
of the full set of n-point Green functions, one is able to extract information on how the 
quantum system compensates for the dissipation due to collisions. We have focused on the 
n-point functions which are likely to survive in the classical limit, i.e. the matrix elements of 
anticommutators of the gauge field operators. Their connected part should correspond to the 
cumulants of the correlation functions in statistical physics. This is one way to characterize 
the fluctuations that maintain the system in thermal equilibrium. 

The result that has been obtained is that the whole series of cumulants differs from zero in 
ultrasoft (or soft) QCD. This result suggests that there may be alternative ways to charac- 
terize these large-scale, classical fluctuations of the (non-abelian) gauge field. 
The third order cumulant has been given special attention. Explicit expressions have been 
written down for the ultrasoft effective theory {p -C gT) and for the soft theory [p <^ T). 
Comparison has been made with the semiclassical approach to the effective theory, where 
a random noise compensates for the dissipation. The solution to the classical equation has 
been written as an expansion in powers of the noise, and the noise average of the 3-point 
function has been computed to lowest (non-zero) order. It agrees with the quantum re- 
sult, which shows that the correspondence is correct between the quantum quantity and the 
classical one. Those results and comparisons have been made at the scale {g'^T In l/g)~^, 
assuming Inl/g ^ 1, i.e. at an intermediate scale between {gT)~^ and {g'^T)~^. This scale 
allows an expansion, in powers of the gauge field, of the Green function {v.D + iC)~^ so that 
explicit expressions have been obtained and compared. 

An extensive use has been made of the special features of the Keldysh basis. The prop- 
erties of the ultrasoft amplitudes may be summarized in the following term in an effective 
action 

r^^\A,,A2) = TT[ [ Ai^{X)j,{X;A,{X))+ [ ]- Ai^iX) }C,.{X,Y; A,{X)) A'^iY) ] (85) 

with = /-c» ^'^ I ^3^- ^i(^) ^2{X) are the two components of the gauge field in the 
Keldysh basis, j^{X; A) and lC^y{X, Y; A) are non-local quantities expressed in terms of the 
propagator (along jagged paths) in the background non-abelian field A. Their expressions 
may be found in Eqs.(PBD(PlD. /C^,^ would vanish if there were no dissipation. Eq.(^) 
explicitly shows that the only existing IPI ultrasoft vertices have either one or two Keldysh 
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indices 2. The Keldysh basis is such that 



A(X) = (A+(X) + A_(X))/2 , A2(X) = A+(X)-A_(X) (86) 

where is the field that fives on tfie forward brancfi of tfie time contour and A_ on tfie 
backward brancfi (in tfie patfi integral formalism). Tfien Ai{X) is tfie average {A^ + A^)/2 
and plays tfie role of a background field, wfiereas A2{X) is tfie difference A+ — A^ and plays 
tfie role of a fluctuating fleld. 

Tfiis interpretation is made stronger if one considers tfie gauge transformation 

A^" ^hA^ - (i/g) h (87) 

In tfiis transformation, tfie propagator in tfie background field A transforms covariantly, as do 
j^(X; A) and /C^,(X, Y; A) , ^ hj^h^ , /C^, ^ hJC^^h^ (D^^j^ = = D^/C^, = D^/C^,) 
Witfi tfie identification A^(X) = A'^(X), rP°\Ai,A2) , defined in Eq.(|8|), wiU be 
invariant under tfie transformation (^) if A2 transforms covariantly 

Ai^^hAi^ h) (88) 



Tfie set of transformation rules (|87D(|88D are tfiose of tfie background field metfiod. Here, tfie 
gauge field is split into a background field ^/^(X) (to be identified witfi tfie average field) 
and a fluctuating fleld a^(X) wfiicfi transforms covariantly under tfie gauge transformation 
of tfie background field, = A^ + a^. One concludes tfiat tfie form (p5D is in tfie class 
of background field gauges 0|. It is wortfi stressing tfiat tfie explicit form of tfie collision 
operator C(v,v') is irrelevant. Its required properties are: symmetry in v v', positive 
eigenvalues, tfie zero mode /^C(v, v') = 0, and gauge-fixing independence. 



Witfi tfie interpretation of tfie A2{X) field as a fluctuating fleld, one is led to tfie idea 
tfiat tfiese fluctuations are, very likely, tfie ones tfiat compensate for tfie dissipation due to 
collisions and maintain tfie system in tfiermal equilibrium. An encouraging fact is tfiat, for 
a classical field, 

- ^ iyFa'iX) }Cf^{X,Y;A) F^\Y) = Sym.Part [j^F'/{X) j;(X; A) = /^j.E(X)) 

(89) 

wfiere Sym.Part means tfie symmetric part in tfie excfiange E\{X) ^ EliY) , and j.E(X) 
is an energy density rate. 
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A Constraint on n-point functions 

and fluctuation-dissipation theorem 

In statistical physics, those theorems relate the fluctuations, i.e. a correlation function's 
cumulant, to the response of the system to an external weak perturbation, and, mostly, they 
involve 2-point functions. In field theory, the response functions are the n-point retarded 
amplitudes (an example is the current of Sec.|) and the statistical correlations may be 
identified with G22-2 at high temperature (see Sec.p.3|). 

There exists relations among the 2" n-point amplitudes which are best stated as a complex 
conjugate relation in the R/A basis. In addition to Trji...r = = Vaa---a , one has, for a 
boson system, [|13|, |10| 



Ar(Pi, . . . Pfc) T{P,A ■ ■ ■ PkA, Pik+1)R . . . Pnn) = 
{-irM{Pu+l ...Pn) r*(Pi, . . . Pfcr, P(,k+1)A ■ ■ ■ Pua) (90) 

where 

k k 

Af{P„ ...P,) = mN{P^) + 1/2) - mN{P^) - 1/2) (91) 

i=l 1=1 

and N{Pi) = (1/2) coth/3p°/2. In particular, A/'(P) = l.(The (-1)" factor depends on con- 
ventions) . 

For the 2-point function, this relation is Tra = ^*ar aiid, with Tjui = 0, the relation 
between the R/A and Keldysh bases written in Sec.^ i.e. Eg. (p!0|) , leads to 

K22{P) = N{P)[Tra - Tar] = N{P) [ K21 - K;, ] (92) 

and to a similar one between the correlation function G22 and the response function Gr 

G'i,{P)=N{P) [Gr-Ga]{P) (93) 

This is a well known fluctuation-dissipation theorem. 

For the case of the 3-point functions, there exists an equivalent relation. Indeed, Eq. (ppp 
relates Trra to F^^^. The relation stemming from Trrji = is in 



- K222 = N{P,)K,22 + N{P2)K2U + N{P,)K221 

+N{P2)N{Ps)K2ii + N{P^)N{Ps)Ki2i + N{P,)N{P2)Kn2 (94) 

Now, Eq. (|To|) and Eq.(^) allow to express K222 in terms of the retarded functions and their 
complex conjugate. The result is 

K222 = N{P2)N{P,){K2u - ^211) + N{P,)N{P3){K,2i - i^i2i) 

+iV(Pi)iV(P2)(Kn2 - K;,,) - (K*ii + ^1*21 + ^n2)/4 (95) 

where the following relation holds for Pi + P2 + P3 = 

N{Pi)N{P2) + N{Pi)N{Ps) + N{P2)N{P3) + 1/4 = (96) 
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K222 is completely determined by the retarded 3-point functions, just as K22 is in terms 
of the 2-point one. However, K22 vanishes if there is no dissipation, while -^^222 does not 
necessarily do so. For a bare vertex (i.e. no dissipation) K222 7^ 0. Indeed, in this case, the 
three retarded functions are identical, such that -ft'211 = ^-^2115 so that relation ( pSf ) becomes 



K'^;^ = <r [-i X 2 + i X 3] = ^K'Z" (97) 

in agreement with Ref.|13[ (where, for a bare n-point vertex, the following result is ob- 
tained: all bare vertices with an even number of 2 vanish, those with an odd number are 
proportional). However, at high T, this bare vertex -ft'222 is lacking a factor T/p^ and is 
subleading compared to the effective vertices. These vertices obey relation (^) with the 
approximation N(Pi) ^ T/p'^ ( and relation ( p6| ) turns into -^^^ + + = 0). It re- 
mains to be seen whether this approximated form of Eq.(^) may be called, in general, a 
fluctuation-dissipation theorem, at high T, for the 3-point function. 

It is worth noticing that, by the argument that leads from Eq.(|^) to Eq.(^), one deduces 
that the correlation function G222 obeys a relation similar to Eq. (^) with the following sub- 
stitutions: G211 = Gr{Pl)Ga{P2)Ga{P3) -^211 Substitutcd tO K211 , 

G2U = Ga{Pi)Gr{P2)Gr{P2) -ft'lii Substituted to K^^^^ , and so forth for the other terms of 



Eq . (p5D . This relation for G222 ^^s obtained in [1^ from the examination of the 3-point 
spectral functions. 

Turning to the 4-point case, F/jk^^ is expressed in terms of F^^^^ij . However, the relation 
between Tjiraa and is only a constraint, and does not determine those amplitudes. 

Consequently, -ft'2222 cannot be expressed only in terms of the response functions, i.e. the 
functions with one Keldysh index 2. More generally, for an n-point amplitude in the R/A 
basis, there are 2^^ — 2 = 2(2^^"^ — 1) non- vanishing amplitudes, with 2"~^ — 1 complex 
conjugate relations, and only n response functions. 

The n-gluon ultrasoft amplitudes are less numerous, because of Bose and gauge symmetries. 
There is one generating functional for the n response functions (one Keldysh index 2), and 
another one for the n(n — 1)/2 functions related to dissipation (two Keldysh indices 2). In the 
correlation functions (which describe the fluctuations) there is an interplay between the two 
sets of functions. The appearance of the same Green function in both generating functionals, 
may be seen as the ultimate fluctuation-dissipation constraint. 

Coming back to the general case, one may tentatively identify the complex-conjugate relation 



(90) with the general constraint of the fluctuation-dissipation theorem. 



B Noise average of classical solutions 

This appendix details the computation of the noise average of the classical 2-point and 
3-point functions examined in Sec.^ One goes to Fourier space. For the noise 

UX,^) = j T&.^"''-'' UP,v) (98) 
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and Eq.(^) translate into 

(( UPu^i) UP2,V2) )) = (27r)V(Pi + P2) Sab ^ 2Civ,,V2) (99) 
For the lowest order solution, Eg . ([7l|) translates into 

= / ^ ^r''{Pl) < ^ V,,\V.P, + tC)-^,^ UPUV2) (100) 

with the convention that a repeated v variable is summed over (compare with Eg . (|37[) ) . 
One now details the computation of the 2-point function. From Eg.( |100|) 

(( A(^Hx)A^'''{y) )) = (^)7 ^-^^''^•^^''^•^H2vr)V(P, + P,) 

where Eg.(p9D has been used. C{v,v') is symmetric in the exchange v <-> v', and so is 
{v.P + iC);;^, . Then, with 2C = [ v.Pi + iC + {-v.Pi + iC) ] /z, one writes 

(i;.Pi + iC)-^ 2C {-v.Pi + tCy^ = [ {v.Pi + iCy^ + {-V.Pi + iC)-^]/i (102) 

Comparing with Eg. fpSD , one recognizes a factor K22 ^v'(Pi) in Eg.(|lOl|). Moreover, smce 

Gr{P2)=Gr{-Pl)=Ga{Pl), 

Gr'{Pi) K22 MPi) ^^''(Pi) = ^ cr' (n. - iia),'u' g:;" = Gr(Pi) (los) 

Pi 

from Eg . (|39|) , so that Eg.( [101| ) turns into 

(( A^)^{X)Af>{Y) )) = / ^ e-^^^ (^-^) 5a, ^ G^^Pi) (104) 



One now turns to the 3-point correlation function. The contribution coming from the 
bare 3-point vertex is obtained from Eg.([79|) 

(( A^-)'^(X)Ar'^(F)4°)''(Z) )) = -Gf{X,Xi) g 

{ {{ Af^''{Xi)Af^^{Y) )) (( (9,,AW^ -9.a(^)(X0AW''(^) )) 

+ {Y,u,b^ Z,p,c) } (105) 



With Eg. (|104|) , one immediately obtains Eg.(^) in momentum space. A similar treatment 



applies to the effective vertex contribution. 

For the third contribution, A^^'^^^X) is in Eq.(^), and one needs, in addition to Eg.( p.04|) , 
another noise average. From Eqs. (|7TD (|55D 

(( UX3,V3) A^'^f{Z) ))=T 6g,G'/{Z,X^) I v^,,Gr{X^,X3-v^,v^) 2G{v^,v;) (106) 
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with the same conventions as in Eg. ([7TD . The 3-point correlation function is, then, in mo- 
mentum space and with the definition ([78|) , 



= -or' (Pi) { (P3) ^Gr{-P,) g mj, T if, 



abc 



+ {P2,i^,b^Ps,p,c) } (107) 



With the identity (|TU^ , and with Gc{-P) = G^{P), the form ([in^) turns into 



X vi^,{v.Pi + iC)-^\^ V2u' [ i-v.P3 + tC)-^ + {v.Ps + iCy^ 

+ iP2,u,b^P3,p,c)} (108) 

One recognizes T(l,2,3) defined in Eg. (p|) in the first term of the expression ( |108| ). The 
second term in (|108D contains 

v,^,{v.Pi + tC)~^l^ {v.Ps + iC)-^l^ v^p, = S{1, 2, 3) (109) 

symmetric in the exchange {Pi, fi' ^ P^^p'). Now, in the symmetric combination ([77|) , 
such terms add up to zero, since 5(1, 2, 3) fabc is cancelled by S'(3, 2, 1) fcba , in contrast to 
r(l,2,3) which obeys Eg.(|4|). 
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